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1. Introduction 

We shall be concerned with the Euclidean correlators of gluons in QCD with an SU (N) local 
gauge symmetry that are fixed to the minimal Landau gauge. These are the fundamental quantities 
in quantum field theory. 

The minimal Landau gauge is obtained by minimizing the Hilbert square norm 



to some local minimum (in general not an absolute minimum) with respect to local gauge transfor- 
mations g(x). These act according to 8 A^ = g~ l A^g + g~ l d^g. At a local minimum, the functional 
F^ig) = \ \ 8 A\\ 2 is stationary and its second variation is positive. It is well known that these two 
properties imply respectively that the Landau gauge (trans versality) condition is satisfied, d - A = 0, 
and that the Faddeev-Popov operator is positive i. e. (©,M(A)©) > for all ft). Here the Faddeev- 
Popov operator acts according to M ac (A)(O c = — d^D a ^{A)(O c , the gauge co variant derivative is 
defined by D a ^(A)(O c = d^co" + f abc A b l co c , and the coupling constant has been absorbed into A. 
Configurations A that satisfy these two conditions are said to be in the (first) Gribov region [[I]] 
which we designate by Q.. It is known that in general there are more than one local minimum of 
F&(g), and we do not specify which local minimum is achieved. This gauge is realized numerically 
by minimizing a lattice analog of F&{g) by some algorithm, and the local minimum achieved is in 
general algorithm dependent. 

The analytic bounds which we shall obtain follow from the restriction of the gauge-fixed con- 
figurations to the Gribov region H, and are the same whether the gluons are coupled to quarks as 
in full QCD, or not, as in pure gluodynamics. In fact the same bounds hold for other gauge bosons 
with SU (N) gauge symmetry, for example, in the Higgs sector, provided only that the gauge-fixing 
is done to the minimal Landau gauge. The numerical results we shall present will be for pure 
gluodynamics in SU (2) gauge theory. 

2. General bounds on free energy 

In the minimal Landau gauge, the free energy W(J) is defined by 



Here the quark degrees of freedom (if present) are integrated out. The Euclidean probability p(A) 
includes the Yang-Mills action, the gauge-fixing factor 8(d - A), the Faddeev-Popov determinant, 
and possibly the quark determinant. We shall use only the properties p(A) > and / dA p(A) = 1. 
The source term J" (x) is real and is taken to be transverse d J = without loss of generality 
because A is identically transverse. The free energy per unit Euclidean volume, w(J) = W(J) /V, is 
the generating functional of connected correlators, 




(1.1) 




(2.1) 



(A(x)A(y)...) 



conn — 



(2.2) 
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The general bound is immediate. From the inequality (J, A) < maxAeci(J,A) = (J, An 
where A max is that configuration in £2 that maximizes (J, A) for fixed J, we obtain 



which gives the bound 



expW(7) < / dAp(A) exp(7,A max ) 
Jo. 

= exp(/,A max ), (2.3) 



W(J) <max(7,A). (2.4) 

AeQ. 



Because Q. is bounded in every direction this bound is finite. It is not hard to show that the 
maximum occurs when A lies on the boundary d£l of CI, and we have the more precise bound on 
the free energy, 

W(J)< max (J, A). (2.5) 

AedQ. 

The right hand side is linear in /, 

max (hJ,A) = h max (J, A) (2.6) 

AedQ. AedQ 

for h > 0. This linear bound is a characteristic feature for integration over any bounded region, 
such as the Gribov region. By contrast, for a free field the free energy is quadratic in J, 

W fKC (J) = (l/2)(J,K- l J), (2.7) 



where K = —d +m . This strongly violates the linear bound (2.5) at large /. 



3. Bound for a plane-wave source 

We now specialize to a plane wave source, 

J£(x) =hcos(k-x 2 )8 a3 S fll , (3.1) 

Here h is the analog in a spin theory of an external magnetic field, modulated by a plane wave. The 
wave number takes on the values k = 2nn/L, where n is an integer, and L is the edge of a periodic 
Euclidean box. The indices 1 and 2 are chosen so J is transverse, d^J^ = 0. For this source, we 
parametrize the free energy per unit Euclidean volume V = L d by 

w(k,h)=W(hcos(k-x))/V. (3.2) 

The gluon propagator is its second derivative at h = 0, 



D(k) = (1/2) 



d 2 w(k,h) 



dh 2 



(3.3) 

h=0 



For a static source, k = 0, we have A max (x) = const, and the Faddeev-Popov operator may be 
diagonalized by Fourier transform. In this case we obtain on a finite volume V = L d , the bound 

w(0,h)<\h\(2n/L). (3.4) 
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Thus on an infinite volume, L — > °°, the free energy vanishes, 

lim w(0,h) = 0, (3.5) 

for a static source of strength h no matter how strong. Lest it be thought that this is a peculiarity of 
the mode k = 0, we also exhibit a bound at infinite volume, with k finite. By an explicit calculation 
one obtains [||] 

w(k,h)<\hk\. (3.6) 



The free energy vanishes in the static limit 



limw(k,h)=0 (3.7) 

k->0 



for all h no matter how strong, in agreement with ( |3.5 



We also obtain a bound on the "magnetization" m(k,h) = Sw jfc ■ Indeed because p(A) is 



normalized, we have w(k,0) = 0, and so 



w(k,h)= / dh'm(k,ti). (3.8) 



o 



Inserting this into the bound (3.6), we obtain, for k > and h > 0, 

1 rh 
h 

which gives 



o 



dti m(k,ti) <k, (3.9) 



l r 

lim - / dti m(k,ti) = 0. (3.10) 

i->0 rt Jo 

Note that m(k,0) = 0, and that dm ^ ^ = - > is positive because w(k,h) is a convex 

function of h. Consequently m(k,h) is positive, m(k,ti) > 0, and we conclude from the positivity 
of the integrand and the vanishing of the last integral that the magnetization produced by a static 
source vanishes, 

limm(k,h) = 0, (3.11) 

for almost all h. Thus, in the minimal Landau gauge, the static static color degree of freedom 
cannot be excited by applying an external color-magnetic field h, no matter how strong. 

Finally we note from the vanishing of w(0,h) = limi ( ^ w(k, h) = that, IF w(0,h) were 
analytic in h, then all derivatives of w(0,h) with respect to h would vanish. In particular the second 
derivative would vanish, - ^ \ h=o = 2D(0) = 0. This is the statement that the gluon propagator 
D{k) vanishes at k = 0. Conversely, if D(0) > 0, then w(0,/i) is non-analytic in h. Numerical data 
indicate that D(0) > in Euclidean dimension d = 3,4, [Q § § [7|, |, g], while D(0) = for d = 2 
JTO] , [I], ll]. 1 Thus non-analyticity of w(0,h) is implied by the lattice data in Euclidean dimension 
d = 3A. 



'The vanishing of D(0) in Euclidean dimension d = 2 is proven in |p^], 
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4. Numerical results for the free energy 

A numerical measurement using a Monte Carlo approach of a current-dependent free energy 
( fOl ) with a current in a fixed gauge is non-trivial, because (yet) no efficient algorithm for Monte 



Carlo updates in a fixed gauge is known, though first proposals exist Q13Q. To circumvent this 
problem, here reweighting will be used. In this case, the lattice configurations are obtained at zero 
external current. The current-dependent free-energy W is then measured by 

e xv(W[J]) = (exp(jjA\\, (4.1) 



which will be performed in the minimal Landau gauge, see [|14J for a review. Since the measured 
quantity is an exponential in the fields, the standard importance sampling cannot be expected to be 
accurate for large J. It will fail at the latest when the source term becomes comparable with the 
action itself. Of course, it cannot be excluded that already a small value of J upsets the importance 
sampling significantly. 



The situation can be a bit better estimated when using the fixed source term (3.1). Using that 



the gauge field in lattice units is bounded by one, the source term is of maximum size 

max J dxJA = Vh (4.2) 
The maximum expectation value of the Wilson action used here is given by 

(43, 

where d is the number of space-time dimensions, = 4/g 2 , and (P) is the plaquette expectation 
value, i. e. the free energy per unit volume. Thus, the maximum h possible is 

6 <M^_2M, (4 . 4) 

which still permits rather large h, under the assumption that going close to this limit will not distort 
the results too severely. In the following some preliminary results from such an evaluation will be 



presented. More results will be presented elsewhere [15]. These will be given for two, three, and 
four dimensions, as it has been found that gauge-dependent quantities depend significantly on the 
dimensionality [Mj. 



With the source (p.l|), the functional W becomes a function of the two independent variables 
k and h. Since k is a lattice momentum, it can have only discrete values, while h is a continuous 
variable. For k, at most ten different values have been used, depending on the lattice size. An 
example for W is shown in figure [[} It is immediately visible that W depends polynomial on h at 
large values of h. Since this sets in already some orders of magnitude below the reliability limit 
(p~4|), this is likely a genuine effect. 

To identify the slope, figure || shows W/Wb, where Wo is \f2nN d ~ l h for k = and N d /y/2hk 
otherwise. This quantity shows the linear dependency on both h and k (ET6|) expected at large h, 



and any correction is found to be smaller than 0(ln(h)) already substantially below the limit (4.4) 
However, the bound Wo is not saturated, and a prefactor smaller than one remains. 
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Figure 1: The free energy density W(h,k)/V in lattice units as a function of h for various values of k is 
shown in the middle panel. The right and left panels show magnifications of the regions at large and small 
h, respectively. The value of h is limited by (4.4), indicated by the green dashed line in the top panel. In the 
right panel the blue dashed line is linear in h, while the red dashed line in the left panel is proportional to h 2 . 
The three-dimensional lattice has volume of 72 3 lattice sites at j3 = 3.73, i. e. (14.4 fm) 3 and a = 0.2 fm. 



Free energy in two dimensions I Free energy in three dimensions I Free energy in four dimensions 




10" 10° 10 2 10' 1 10 10 2 10* 10° 10 2 10' 1 10 10 2 10' 10 3 Iff 1 10 1 1 10 10 2 

h h h 



Figure 2: The normalized free energy W/Wo as a function of h for two (left panel, 240 2 /(48 fm) 2 at j3 = 
1.991a = 0.2 fm), three (middle panel, 36 3 /(7.2 fm) 3 at j3 = 3.73/a = 0.2 fm), and four (right panel, 24 4 /(4.8 
fm) 4 at j3 = 2.221/a = 0.2 fm) dimensions. See jl5[ ] for an explanation of the somewhat unusual error 
distribution. The different colors and symbols refer to increasing values of k, with values depending on the 
lattice parameters, but generically the large the value of W/Wo, the smaller k and the largest values always 
correspond to k = 0. 
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Figure 3: The value of the normalized free energy density W(h,k)/Wo at k — in the domain at large h, 
where it becomes constant, as a function of lattice volume and discretization for two (bottom left panel), three 
(bottom right panel) and four (top left panel) dimensions. The top right panel shows a cut at fixed spatial 
volume. Statistical error bars are smaller than the symbols. The volume is (12 fm) 2 in two dimensions, (3.6 
fm) 3 in three dimensions, and (1.2 fm) 4 in four dimensions. 



To investigate whether this is a lattice artifact and/or depends on the dimensionality, the re- 
maining constant of proportionality for a range of lattice sizes and discretizations has been deter- 
mined [15], and preliminary results are shown in figure At first sight, no qualitative difference is 
found between different dimensions. It is visible that at fixed lattice spacing the expected bound is 
less fulfilled the larger the volume. However, when moving towards the continuum at fixed volume 
the ratio stays, more or less, constant. No final conclusion can be drawn from this, except that the 
order of limits could be quite important. 

Concluding, the numerical results suggest that the qualitative behavior appears to be indeed 
the one expected from the analytical analysis of section |3| Given that this behavior sets in already 
several order of magnitudes in h below the limit (|4.4|), there is a fair chance that this is a genuine 
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effect and not an artifact of reweighting. However quantitatively the agreement is less convincing, 
and the bound on the free energy is significantly undersaturated. Since a saturation is only expected 
asymptotically, this may, or may not, be an artifact of reweighting. It may, of course, be that 
configurations saturating the bound are not sufficiently sampled. This may either be due to the 
importance sampling, or because minimal Landau gauge usually misses the Gribov copies closest 



to the horizon []14[]. Some further considerations will be presented elsewhere Q15Q. Furthermore, 
the dependence on the lattice parameters seen in figure |] suggests that the order of limits may be 
important, which could indicate that non-analyticities may play a role. All of this recommends 
further investigations, which are currently under way. 
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